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ABSTRACT 

Keck near-infrared images of Neptune from UT 26 July 2007 show that the cloud feature 
typically observed within a few degrees of Neptune's south pole had split into a pair of bright 
spots. A careful determination of disk center places the cloud centers at —89.07 ± 0.06° and 
—87.84 ± 0.06° planetocentric latitude. If modeled as optically thick, perfectly reflecting layers, 
we find the pair of features to be constrained to the troposphere, at pressures greater than 0.4 
bar. By UT 28 July 2007, images with comparable resolution reveal only a single feature near 
the south pole. The changing morphology of these circumpolar clouds suggests they may form 
in a region of strong convection surrounding a Neptunian south polar vortex. 
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1. Introduction 



The atmosphere of Neptune is active. In ad- 



in the stratosphere and troposphere (Orton and 



dition to global changes on decadal scales (Lock- 



wood 



wood and Jerzykiewicz 2006 Hammel and Lock- 



2007), individual clouds evolve on time 



scales as short as hours (e.g. Limaye and Sro- 


movsky|1991||Sromovsky et al.|1995, 


2001). Near- 



infrared (NIR) images, which probe Neptune's 
lower stratosphere and upper troposphere in re- 
flected sunlight, reveal bright bands and cloud 
features against a dark disk. The pattern of zonal 
circulation, which is most pronounced at mid lati- 
tudes, continues to high southern latitudes, where 
a bright, unresolved cloud feature is present within 
a few degrees of the pole. This feature has been 
seen since the Voyager 2 era dSmith et aL]|1989[ ). 



Interest in Neptune's global dynamics has been 
reignited by mid-IR observations of enhanced 
emission over the south pole relative to the rest 



of the planet (Hammel et al. 2007 Orton et al. 



2007). This enhancement suggests a temperature 
increase of 4-5 K in the south polar region, and 
has led these authors to draw parallels with Sat- 
urn, where the summer pole is similarly warm 
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Yanamandra-Fisher 2005). Recently it has been 



observed that both of Saturn's poles exhibit a 
localized hot spot surrounded by cooler zones, 
forming circulation cells that make up part of a 
planet-wide circulation ( Fletcher et al.|20Q8 Dyu- 



dina et al. 2008). A global dynamical pattern 



has also been inferred for Neptune's atmosphere 



(Martin et al. 2008); it remains to be seen if this 
pattern includes a polar circulation cell like those 
on Saturn. 

We obtained high spatial resolution NIR images 
of Neptune's south pole in July 2007, revealing 
temporal evolution of discrete circumpolar clouds. 
Here we describe these features, model their prop- 
erties, and discuss the implications for the dynam- 
ics of the Neptunian south polar region. 

2. Observations and Data Processing 

We observed Neptune from the 10-meter W.M. 
Keck II telescope on Mauna Kea, Hawaii, on 26-28 
July (UT) 2007, as part of a long-term program 
studying the planet's atmosphere in the NIR. J-, 
H-, and Kp-band images were taken using the nar- 
row camera of the NIRC2 instrument, coupled to 
the adaptive optics (AO) system (Table 1). The 
1024x1024 array has a pixel size of 9.963 ± 0.011 
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mas/pixel in this mode f Pravdo et al.|20Q6 ), which 
at the time of observations corresponded to a phys- 
ical size of ~210 km at disk center. 

The data were flat fielded and sky-subtracted, 
and bad pixels were replaced with the median of 
the 8 surrounding pixels. All images were cor- 
rected for the geometric distortion of the array 
using the 'dewarp' routines provided by P. Brian 
Camerorj^ who estimates residual errors at ^ 0.1 
pixels. We measured a full width at half maxi- 
mum (FWHM) of 0.039 ±0.005'' for a stellar point 
source (SAO 146732) on the days of observation, 
which is consistent with the diffraction limit of the 
telescope at 2 /im, and corresponds to an effective 
resolution of ~800 km at the center of the disk. 

The images were photometrically calibrated us- 
ing the star HD 22686, then converted from units 
of observed fiux density to units of I /F^ which is 
defined as ( Hammel et al.|p^89 ): 



F 



r^F^ 

n Fq 



where r is Neptune's heliocentric distance, ttFq is 
the sun's fiux density at Earth's orbit (Colina et al 
1996), Fn is Neptune's observed fiux density, and 
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Q is the solid angle subtended by a pixel on the 
detector. By this definition, //F = 1 for uniformly 
diffuse scattering from a Lambert surface when 
viewed at normal incidence. 

The 26 July data (Fig. [T]) show two distinct 
features near Neptune's south pole. Other unre- 
solved sources, such as moons and other clouds, 
do not appear double, so it seems unlikely this is 
an artifact. The 28 July data (Fig. |2|, at com- 
parable resolution, show only a single unresolved 
polar cloud feature. To determine the total I/F 
of the unresolved features above the background, 
we model them as one (for the single cloud on 28 
July) or two (for the pair of clouds on 26 July) 2D 
Gaussians. We use the best-fit parameters from 
the least-squares fitting routine MPFITQ to find 
the I/F of each cloud feature (Table 2). From 
experience we estimate the error in the photom- 
etry to conservatively be 20%; the errors in the 
cloud fits are between 4% and 12%. In Kp band, 
we do not see the polar cloud features above the 
noise (Fig. [T]). To set an upper limit for the I/F 
of the clouds in Kp band, we assume the clouds 
will be detectible when the mean value of the sig- 
nal plus noise is equal to one standard deviation 
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Fig. 1. — J-, H- and Kp-band images from 26 July 2007 in units of I/F. The first three images in each filter 
were averaged together to produce the figure. The top panels all have same color bar. In the top images, 
the pole is bright due to the convergence of longitude lines; therefore we zoom in on the pole in the bottom 
panels. The top row highlights the features by superposing contours at 10% 20% and 26% above the H-band 
background; and 3%, 5% and 10% above the J-band background level. Latitude circles at 4° spacings are 
also shown. The arrows indicate the positions of the two cloud features; in Kp band the clouds are not above 
the noise; therefore no intensity contours are plotted. 



above the noise mean. We confirmed this upper 
limit by convolving the maximum Kp I/F with a 
Gaussian having the FWHM of the point spread 
function (PSF) and adding the observed level of 
noise. 

3. Navigation and Cloud Locations 

Neptune's bright and dynamic activity limits 
the accuracy of image centering when using the 
planet's limb as a reference, since the PSFs of 
bright features near the limb often extend off the 
edge of the disk. To perform the centering more 
precisely, we found the positions of three of Nep- 
tune's moons in each H-band image, and shifted 
the images in x and y so that the positions of the 



moons were best fit to lie on their orbits, as de- 
rived from the Planetary Rings Node ephemeris 
dataj^ (Fig. [3|. While other methods of aligning 
images do not necessarily give you the location of 
disk center, this technique allows us to improve 
both the relative alignment of the images as well 
as the determination of disk center coordinates. 
The centering error in our least-squares fit is on 
average 0.4 pixels per image, which agrees well 
with the image-to-image scatter of the positions 
of individual cloud features. This is an improve- 
ment over our centering by limb detection, which 
we estimate to result in 1-2 pixel errors in the de- 
termination of disk center. In J- and K- band we 
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Fig. 2. — Average of first three H-band images 
from 28 July 2007, in units of I/F, with latitude 
and longitude lines plotted. The middle image 
highlights the single cloud feature with contours 
of 5% 10% and 20% above the background level. 



do not have a sufficient number of images to im- 
plement our moon centering technique; therefore 
these images are aligned to the H-band images by 
eye, to sub-pixel accuracy. 

To find the positions of the clouds near the pole, 
we use the moon-centered H-band images and the 
solutions from the 2D Gaussian fits (Section 2). 
We determine the weighted mean location of each 
cloud in image coordinates, which we then trans- 
form into planetocentric latitude and longitude, 
using the JPL Horizons ephemeris informatior]^ 
for sub-observer latitude and longitude at the time 
of the observations. The errors in the latitude 
and longitude of the clouds are estimated using a 
Monte Carlo method, whereby for each parameter 
involved in calculating the latitude and longitude. 
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a random distribution of values consistent with its 
probability distribution is generated. These pa- 



rameters include the array rotation angle (Pravdo 



et al. 2006), image center coordinates, and cloud 
center coordinates. For each set of generated pa- 
rameter values, the latitude and longitude are then 
calculated; errors are determined as the standard 
deviations of these simulated data sets. 

We find that on 26 July the two cloud features 
are near, but not directly at the pole as determined 
by moon orbits: the spot nearest the pole is lo- 
cated at -89.07 ±0.06° latitude and 293 ±3° lon- 
gitude. The second cloud resides approximately 
500 km further from the pole, at -87.84 ± 0.06° 
latitude, 238 ±2° longitude. The single spot on 28 
July is likewise near, but not at, the south pole, at 
-88.47 ± 0.03° latitude, 86 ± 2° longitude (Table 
2). 

4. Radiative Transfer Modeling of Fea- 
tures 

The J, H, and Kp filters probe various depths 
through the stratosphere and upper troposphere 
(Fig. |4|; therefore we can use the data to set lim- 
its on the heights of the south polar cloud fea- 
tures. To do this, we adapted the radiative trans- 
fer (RT) code for Titan from 



Adamkovics et al. 



( [2QQ7| ) for Nept une. We adopt the Neptune tem- 



perature profile derived by Lindal ( 1992 ) from oc- 



cult at ion measurements, and incorporate the ma- 
jor opacity sources that dominate Neptune's 1- 
5 /im spectrum: we use the H2 collision-induced 
absorption coefficients for hydrogen and helium 



(Borysow et al. 1985, 1988 Borysow |1991| |1992 



1993), assuming an equilibrium ortho/para ratio 



for H2; and the k-distribution coefficients for H2- 



mg 



broadened methane ( Irwin et al.^^2006) . Follow- 



Gibbard et al. (2002), our lower boundary is 
an optically thick cloud at 3.8 bar, presumably 
the H2S cloudtop (de Pater et al. 1991), with a 



Henyey-Greenstein asymmetry parameter of —0.1. 
Aerosols are treated as Mie scatterers with real in- 
dices of refraction of 1.43 and imaginary indices 
of refraction = 0. Within a given haze layer, 
aerosols are characterized by a single particle ra- 
dius; this radius, along with the altitude range 
and number density of aerosol particles, are free 
parameters, chosen to fit the data. The RT equa- 
tions are then solved using a two-stream approxi- 
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Table 2: Spot locations and flux densities in units of I/F, as well as the I/F of the region surrounding the 
clouds. For the clouds, the maximum Kp/J and Kp/H ratios allowed by the data are also given. 
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0.008 ±0.002 
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28 July 2007 














-88.47 ±0.03° 


86 ±2° 


0.017 ±0.009 


0.042 ±0.009 


< 0.005 


0.3 ±0.1 


0.12 ±0.01 


background 




0.014 ±0.003 
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mation ( |Toon et al.||1989D . 

To model the south polar cloud features, we 
first developed two cloud-free 'background atmo- 
sphere' models. For the first model, we choose the 
simplest haze distribution consistent with the J-, 
H-, and K-band I/F values in the region surround- 
ing the polar clouds (Table 3). For comparison, we 
also find a best model fit to a medium-resolution 
spectrum from the IRTF Spectral Library (Fig. 
[5|, taken with the SpeX instrument on the NASA 
Infrared Telescope Facility (IRTF) on Mauna Kea 



(Rayner et al. 2009). These data are averaged over 
a significant fraction of Neptune's disk. Although 
this disk-integrated model has little contribution 
from the south pole, this second model, which has 
very different haze properties from the first (Ta- 
ble 3), provides a qualitative test of the sensitiv- 
ity of our cloud modeling to the choice of model 
hazes. In addition, fitting to the SpeX data al- 
lows us to model the albedo of the deep 'surface' 
cloud: we find that given our choices of bottom 
cloud depth and scattering properties, the single- 
scattering albedo {u) of the cloud must decrease 
from 0.55 at 1.1 fim to 0.20 at 1.27 jam and 0.13 at 
1.56 /im. These values are reasonably consistent 



with |Roe et al.| ( |2001 ). 

We then model each of the south polar clouds 
as a perfectly reflecting {uj = 1), optically thick 
cloud layer with a Mie asymmetry parameter of 
0.5. We consider clouds at each of 100 altitudes 
between 0.4 mbar and 3.8 bar. For each case, we 
average the resulting 'cloudy' model over H, J, and 
Kp bands for comparison with our image data, and 
subtract the I/F of the 'cloud- free' model to get 
the I/F contribution from the cloud. 

For the 26 July data, we find the ratio of the 



Table 3: Haze layers, chosen to fit the I/F in the 
three spectral windows in the the region around 
the polar clouds (top); and the IRTF SpeX data 
(bottom). 
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intensities of the two clouds, which is affected 
by spot fitting errors but not calibration errors, 
to be 1.1 ± 0.1 in H band and 1.0 ± 0.1 in J 
band. Therefore, our results are consistent with 
the clouds being of the same size and at the same 
altitude. However, since the clouds are unresolved, 
we do not know the pixel filling fraction for each 
of the clouds. We can set an upper limit on 
the clouds' altitudes by comparing the maximum 
model-derived Kp-to-H and Kp-to-J I/F ratios to 
the observed values: the higher a cloud is in the 
atmosphere, the greater the expected Kp-band in- 
tensity. The upper limit Kp-to-H and Kp-to-J in- 
tensity ratios for each of the clouds, given the un- 
certainties in H and J, are presented in Table 2. 
Figure [6] shows the results of the modeled Kp/H 
and Kp/J intensities, compared to the maximum 
of these ratios from the 26 July data. 

We find that the upper limit to the altitude 
of both clouds on 26 July is 0.4 bar. This also 
provides an estimate on the minimum size of the 
clouds: if either of the clouds is at 0.4 bar, it must 
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Fig. 3. — Fitting disk center using Neptune's 
moons. The positions of Galatea (green), Larissa 
(blue) and Despina (red) were found in each image 
and aligned with the moon orbits to establish the 
planetary coordinates. The top panel shows the 
moons in the first H-band image on 28 July. The 
middle plot shows the orbits relative to the posi- 
tions of the moons in the images when centered 
using the planet's limb. The bottom plot shows 
the moon positions after fitting them to their or- 
bits. 



fill at least 10% of a NIRC2 narrow camera pixel 
to be consistent with the observed H and J-band 
intensities. A deeper cloud would necessarily be 
larger to produce the same values of I/F. For the 
single cloud on 28 July, we find that its altitude 
also has an upper limit of 0.4 bar. 
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Fig. 4. — Two-way transmission in Neptune's 
atmosphere as predicted by our RT model. The 
curves show the pressure levels at which 10% and 
90% of the light is transmitted to the top of the 
atmosphere and back. 
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Fig. 5. — Comparison of radiative transfer model 
(green) and long-slit IRTF SpeX spectrum (black). 



The models used in Fig. [6] use the haze dis- 
tribution that best matches the I/F of the back- 
ground atmosphere near the cloud features. How- 
ever, we expect our choice of haze parameters to 
have little effect on our results. This is because we 
are interested primarily in the I/F increase due to 
adding a cloud into the background 'cloud-free' at- 
mosphere. When we ran our second set of models 
using the simpler haze distribution used to fit the 
IRTF SpeX data (Table 2), we found that chang- 
ing the haze parameters affected the 'cloud-free' 
and 'cloudy' spectra in a similar way, so that to 
within 3%, the calculated change in I/F from a 
given cloud was independent of the haze parame- 
ters. 

5. Discussion 



A small bright feature at Neptune's south pole 
has been consistently observed since the Voyager 
era, implying a stability that contrasts the short 
(^hours) lifetime of most small cloud features 
at lower latitudes ( Limaye and Sromovsky|[l991 ). 
Our derived upper limit of the altitude of the cir- 
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Fig. 6. — These plots show the model results 
for Kp/H intensity ratio (solid red) and Kp/J ra- 
tio (solid cyan) as a function of cloud altitude for 
each of the two cloud features from 26 July 2007. 
Dotted lines are at the maximum ratios allowed 
by the data, showing that the Kp/H ratio sets a 
stronger limit on the maximum cloud height than 
the Kp/J ratio. The allowed parameter space is 
cross-hatched. 



cumpolar features observed here is consistent with 
clouds formed by the upwelling and condensation 
of methane gas, suggesting this may be a region 
of continuous cloud formation by convection. 

The persistence of south polar cloud activity 
suggests that there is an organized circulation pat- 
tern at Neptune's south pole. Such a pattern has 



been inferred before: Hammel et al. (2007) show 



that there is a mid-IR temperature enhancement 



at Neptune's south pole, and Martin et al. (2008) 



find that the pole is warmer than its surroundings 
at microwave wavelengths as well. The latter au- 
thors suggest that a global dynamical pattern, in 
which air near the south pole is subsiding in both 
the stratosphere and troposphere, can explain the 
observed high temperatures at both the mid-IR 
and radio wavelengths. Subsiding motions will 
adiabatically heat the atmosphere, visible as el- 
evated temperatures in the mid-IR. The subsiding 
air is likely dry, as condensable gases will condense 



out during the ascending branch of the dynamical 
pattern. Such dry conditions enables one to probe 
much deeper warmer levels at radio wavelengths, 
leading to temperature enhancement in microwave 
observations. 




Fig. 7. — Map of Saturn's SPV [adapted from 



Dyudina et al. (2009)] illustrating the small, bright 
clouds encircnng the storm's 'eye', many of which 
have the same general appearance between 2004 
and 2006 observations. Dyudina et al. (2008) con- 
clude these clouds are located in ISaturn's lower 
troposphere. 

On Saturn, a south polar circulation cell has 
been observed, in the form of subsidence in a po- 
lar vortex surrounded by a region of upwelling 



( [Fletcher et al.||2008[ ). The structure of Saturn's 
south polar vortex (SPV) ( [Dyudina et al.||2009l ) 
possesses similarities with terrestrial hurricanes, 
such as a well- formed central eye, concentric eye- 
walls, and a surrounding ring of strong convection 
(Fig. [t]). The region corresponding to the eye of 
the SPV is a hot, nearly circular region within a 
diameter of 4200 km around the south pole. The 
eye is mostly clear of clouds above ~1 bar (Dyu- 



dina et al. 2008), and depleted in phosphine gas 
( Fletcher et al.[[20Q8 ), suggesting that this is a re- 
gion of subsiding air. Surrounding the eye is a ring 
of discrete, bright clouds in the lower troposphere, 
which may be analogous to the heavily precipitat- 
ing clouds encircling the eye of terrestrial hurri- 
canes (Dyudina et al. 2009). The similarity be- 



tween the temperature enhancement at Neptune's 
south pole to the enhancement at Saturn's south 



pole (Hammel et al.[[2007, and references therein) 
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first led to the suggestion tiiat Neptune, like Sat- 
urn, harbors a long-lived hot south polar vortex 
(Orton et al. 2QQ7| ). Our finding that Neptune's 
south polar cloud feature (s) may be indicative of 
rigorous convection near, but not at, the pole, sup- 
ports the analogy between Neptune and Saturn's 
south polar environments. 

The unexpected discovery of a vortex at Sat- 
urn's north (winter) pole ( Fletcher et al.|2008 ) im- 
plies that polar vortices can exist despite consid- 
erable variations in seasonal insolation, and may 
be general features of giant planet atmospheres. 
However, possible analogues to Saturn's vortices 
have yet to be directly observed, as they await 
the measurement of polar wind speeds by high- 
inclination space missions. Further details about 
the morphology of clouds near Saturn's north pole 
may help us understand the types of cloud mor- 
phologies we might expect around polar vortices, 
and provide insight into how to interpret the be- 
havior of Neptune's polar cloud features in the 
context of a possible vortex. 

6. Summary 

We have observed a transient double cloud fea- 
ture near Neptune's south pole on 26 July 2007. 
The locations of the two features, as determined 
by the orbits of three of Neptune's moons, suggest 
they are both circumpolar clouds separated in lat- 
itude by ^ 1.2°. The single circumpolar feature 
seen on 28 July 2007 is at an intermediate latitude. 
Radiative transfer modeling indicates that the fea- 
tures on both days are at depths of greater than 
0.4 bar, which is consistent with the formation of 
methane condensation clouds. The morphological 
change we observe on 26 July, combined with the 
location of the single cloud on 28 July, may indi- 
cate that the bright south polar spot seen since 
Voyager is not a single stable feature, but rather 
the signature of persistent cloud activity related to 
strong convection, perhaps in the region surround- 
ing a polar vortex. Continued high-resolution ob- 
servations of Neptune's south pole will provide 
more insight into the dynamics of this region. 
Such observations should be performed over sev- 
eral hours to capture the evolution of the features. 
We also await more high-inclination observations 
of planetary atmospheres to provide further con- 
text for understanding giant planet polar environ- 



ments. 
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